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ABSTRACT 

We consider the problem of an elliptical background source crossing a linear gravitational lensing 
fold caustic. We derive a simple expression for the light curve of a source with a uniform surface 
brightness that is accurate to third order in the ellipticity e (yielding brightness errors of < 1% for 
e < 0.3). We then consider caustic crossings of a rotating star, an oblate giant extrasolar planet, and 
an inclined standard thin accretion disk around a quasar black hole, and find the following results. (1) 
If most stars have rotation periods similar to chromospherically-active (i.e. spotted) giant stars in the 
bulge, than ~ 15% should be sufficiently oblate (e > 0.25) to produce detectable (> 1%) deviations 
with current observations. The form of the deviation due to ellipticity is qualitatively similar to 
that due to limb-darkening. Thus stellar oblateness, in general, may have to be taken into account 
in interpretations of precise limb-darkening measurements with microlensing. (2) Giant planets will 
generally not produce a detectable oblateness signal, either because they rotate too slowly (if they 
are close-in and tidally locked to their parent star), or because they are too faint for detection (if 
they are farther away from the parent star). However, a close-in planet with a large ellipticity could 
cause a detectable distortion, especially if the caustic crossing occurs in the direction parallel to the 
major or minor axis of the planet's surface. (3) There is a near-degeneracy between ellipticity and 
position angle for equal-area sources with scale-free intensity profiles and elliptical isophotes. This 
degeneracy results in a factor of '--^ 2 uncertainty in the measurement of the scale-length of a standard 
thin accretion disk using observations of a single, linear fold caustic-crossing quasar microlensing event 
in passbands probing the outer parts of the disk. More precise or higher-frequency observations can 
reduce this uncertainty. 

Subject headings: gravitational lensing - stars:atmospheres - planetary systems 



1. INTRODUCTION 

Gravitational lens caustics have proven enormously 
useful in studies of a broad range of astrophysical phe- 
nomena. Caustics are closed curves that describe the 
set of source positions where the Jacobian of the map- 
ping from source plane to image plane induced by the 
lens vanishes. The magnification of a point source on 
a caustic in formally infinite. Caustics are composed 
of two different types of singularities: folds and cusps. 
Folds are smooth, curved lines that meet at cusp points. 
Cusps are higher-order singularities where the tangent 
to the caustic curve is undefined. The majority of the 
length of a typical caustic curve is well-described by a 
simple fold singularity, and therefore fold caustics are 
generally more common. A source sufficiently close to a 
fold caustic is lensed into a pair of equal-magnification, 
opposite-parity images whose magnification diverges as 
(Am_l)~^/^, where Au_l is the perpendicular separation 
of the source from the caustic. These are the two images 
that are created locally by the fold caustic. There will 
generally be other images created by the global poten- 
tial of the lens that are not associated with the caustic. 
However, these images will typically behave in a smooth 
and continuous manner near the the caustic. Here we 
consider the case of microlensing of sources near fold 
caustics. In this case, the individual images are unre- 
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solved, and only the total magnification of all the images 
is observable. Generally, microlensed sources move on 
relatively short time scales, and therefore the change in 
magnification as a function of time, i.e. a light curve, is 
observable. Due to the divergence in magnification as a 
source crosses a fold caustic, the light curve can be used 
to study sources at higher spatial and angular resolution, 
than would be possible with conventional techniques. 

Caustics have been proposed or implemented to 
study objects ranging in scale from planets to quasars. 
Applications relevant to the discussion here in- 
clude the direct detection and characterization of 
close-in extrasolar giant planets and associated struc- 
tures via their reflecte d ligh t (Ashton & Lewis 200 i| 
Gaudi. Cha ng, fc Hani l2003t iGraff'fc Gaudil 2000; 



Lewis &: Ibata.20p0ll. the precise measurement of stellar 
limb- d arkening llAbe et all 120031 lAlbrow et all 119991 
120011 [Fields et al.' '200^, and the resolution of the 
central engines of quasars ( e.g. lA gol & Krolik 199^ 
Goicoechea et all I20'03t iGo uld & Miralda-EscucQ 
1997t iGrieger. Kavser & Refsdal 1981 



-sfiT- A/; Ke 

Grieger. Kavser. fc Schramml Il991j) . We note that, 
for systems composed of point masses, extended 
caustics occur only for composite lensing bodies, the 
simplest example being a binary star. A single star 
also produces formally divergent magnification, at 
the point at which the observer, lens, and source are 
perfectly aligned. Such fortuitously aligned systems 
can then be used to resolve stellar surfaces. Indeed, 
microlensing by a point mass h as been proposed as a 
— ^— i tool to study limb darkening (iGaudi fc Gouldl 119991: 
^IHendrv et alT iHevrovskvlWo,^ l\ 
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1998), stellar spots llllendrv. Brvce. fc Valls-GabaudI 
2002t iHevrovskv fc Sasselovl \20m . and ellipticity 
ijHevrovskv fc Loeblll997|) . The lens and source must be 
aligned to better than the angular size of the star for the 
source to be resolved in single-lens events. However, the 
source size for typical Galactic microlensing events is less 
than a few percent of the angular Einstein ring radius of 
the lens. Therefore, caustic-crossing binary-microlensing 
events are more common than source transit events by 
single stars. 

Analytic results for the magnification near fold caustics 
have been derived only for circularly-symmetric sources, 
and only for a few special functional forms for the surface 
brightnes s profile of such sources (see. e.g.fAlbro\v et al.l 
IT999t iDoi^ihlig 1200^ ISchneider. Ehlers fc Falcol 119921) . 
However, for a number of applications, one is interested 
in more complex source geometries. Here we consider 
one such geometry that should occur frequently in na- 
ture: elliptical sources. In fjJlwe derive a simple, analytic 
expression for the change in magnification relative to a 
circular source that is accurate to third order in the ellip- 
ticity. We apply our results in fj^J considering the mea- 
surement of stellar oblateness and limb-darkening using 
fold caustic crossings f ij3.1|l . the possibility of measuring 
the oblateness of giant extrasolar planets ft |3.2|l . and the 
resolution of a quasar accretion disk We summa- 

rize and conclude in ^ 

2. ELLIPTICAL SOURCES NEAR FOLD CAUSTICS 

The total magnification of a point source with an an- 
gular sepa ration of Au \ normal to a simp le linear fold 
caustic is (jSchneider. Ehlers fc FalcdlT99^ . 



Au I 



^1/2 



H{Auj_), 



(1) 



where Ur is related to the der ivatives of the 
lens mapp ing at the fold fsee iGaudi fc PettersI 

2onl 



It 



Fetters. Levine, fc \ yambsgansl l2001t 
Schneider. Ehlers fc Falco 199^^ and H{x) is the 



Heaviside step function. Here and throughout, all angu- 
lar separations on the sky will be in units of the angular 
Einstein ring radius of the lens, 9e — (2i?s/-D)^/^, 
where i?s — 2GM/c^ is the Schwarzschild radius of mass 
M, D — DaiDog/ Dis, and I?os, D^s, and D\s are the 
distances between the observer-lens, observer-source, 
and lens-source. For cosmological objects, angular 
diameter distances should be used. 

The magnification of a finite source near a fold is just 
the convolution of equation (^) over the source. 



A,. 



/^du/(u)A(u) 
^du/(u) 



(2) 



where /(u) is the surface brightness distribution of the 
source, and D is the solid angle extended by the source. 

For the simple case of a uniform-surface-brightness 
source, adopting a polar coordinate system with origin 
at the center of the source (see Figure^, equation |j2Jl 
can be rewritten as 

1/2 -, i-2n i-r{e;e.tl>)/p 

d9 dx 
'o Jo 



[t] + xcos 6']i/2 



(3) 




Fig. 1. — The geometry of the integration of the magnification 
of an eUiptical source near a fold caustic in a polar (r, 9) coordinate 
system. The origin is chosen to be the center of the ellipse. The 
major axis of the ellipse is indicated by the dashed line, which is 
at an angle 4> with respect to the perpendicular to the caustic. 



Here, p is the characteristic angular size of the source, 
defined such that the total area extended by the source 
is irp^ (i.e. p — r for a circle, and p — \fab for an 
ellipse); r(9;e,(f)) is the outer boundary of the source. 
We have also defined the dimensionless angle x = u/p, 
and distance from the caustic, rj = Au±/p. The inner 
integral in equation ^ can be evaluated analytically, 
although care must be taken to establish the limits of 
integration so that the condition rj + xcos 9 > is met. 

For a circular source, r = jO, and the magn ificati on has 
the w ell-known form (see, e.g. ISchneider. Ehlers fc Falcol 
IT991 . 

^('7,0) -^-j Goiv), (4) 



where 



Gniv) 



^1/2 



(n + 1)! 
{n + 1/2)1 



dx 



max( — 77, — 1) 



(1 



^2-)n+l/2 

{x + 77)1/2 



(5) 



The function Go{ri) is shown in Figure ^ 

We now consider a uniform source with arbitrary ellip- 
ticity e. We first address the case of r? > 1, i.e. prior to 
the source's entry onto the caustic, for which the limits 
of integration of the inner integral in equation Q are 
(0,r/p) and the limits of the outer integral are (0, 27r). 
The magnification is then. 



P 



1/2 



27r 



dO- 



x{[(r//9) cos6' — 2ri\[{r/ p) cos ( 



3 cos2 9 
77]i/2 + 2773/2}. 



(6) 



For an ellipse with semi-major axis a, r = a[cos a -I- 
(1 — e^)^"'^ sin^ aj^-'^/^, where a = 6 — (j), and (j) is the 
angle of the major axis with respect to the normal to the 
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Fig. 2. — (a) The dotted line shows the function the 
analytical approximation describing the normahzed difference in 
magnification between an equal-area elliptical and circular uniform 
source as a function of the angular separation r\ between the cen- 
ter of the source and the caustic in units of the angular size of 
the source p. For sources of ellipticity e and position angle (/i, this 
curves should be multiplied by e'^cos{20). See equation ISl The 
dashed line shows Gxj^^ii) — Goi'']), which is the function describ- 
ing the difference in magnification between a uniform and limb- 
darkened circular source, (b) The solid line shows the normalized 
light curve Go{ri) for a uniform, circular source crossing a fold 
caustic. The dotted line shows the exact light curve for a uniform 
elliptical source with ellipticity e = 0.5 and (f> = 90°. The dashed 
line shows Go + r[Gi/2{*?) — Go{v)]y light curve for a circular, 
limb-darkened source with limb-darkening parameter F = 0.5. For 
all curves in both panels, the magnification for a source of size p 
can be found by multiplying by p^^^^. 



Fig. 3. — In all panels, the dotted line shows the function 
cos(2<f>)E{r^), the analytical approximation describing the differ- 
ence in magnification between an elliptical source with ellipticity e 
and position angle ij>, and an equal-area, circular, uniform source, 
as a function of the angular separation rj between the center of the 
source and the caustic in units of the angular size of the source 
p. The solid line shows the exact solution for the same quan- 
tity, obtained numerically from equation |^ The dashed line shows 
r[Gi/2('?) — Go{il)]y the difference in magnification between a cir- 
cular, limb-darkened source with limb-darkening parameter F, and 
a circular, uniform source. Here F = — 2e^ cos(2(/)), which approx- 
imately reproduces the deviation due to ellipticity. For all curves, 
the magnification for a source of size p can be found by multiply- 
ing by p-i/2. (a) e = 0.5, <f> = 90°, and F = 0.5. (b) e = 0.3, 
= 90°, and F = 0.18. (c) e = 0.05, </> = 90°, and F = 0.005. For 
ellipticities as small as shown in panel (c), the analytical approx- 
imation becomes near-exact, and the dotted and solid curves are 
indistinguishable. 



caustic (see Figure ^ for an illustration). Note that we 
have centered the coordinate system at the center (rather 
than at the focus) of the ellipse. An ellipse with area np'^ 
has a = p{l — e^)~^/^, and thus we have, 



[cos^ a -f- (1 - 
1 + ie^ cos2a-|-0(e'^), 



-isin^d-i/^ 



(1 



-1/4 



(7) 



Inserting equation (jT)) into equation ©, keeping terms 
only up to order e^, and discarding terms that are odd 
under 9 —0, we find after some algebra, 



Ai7^;e,^)=A{ri.,0) 



1/2 



cos{2(f>)E{r/), (8) 



where we have defined the function 



1 

47r 



d0 



cos 26* 



(cos6i-|-r;)i/2' 



(for 77 > 1) (9) 



which can be expressed in terms of elliptic integrals. 

These results assume that the source is off the caustic. 
During caustic crossing, when \r]\ < 1, the limits of in- 
tegration become complicated. However, in the limit of 
e — *■ 0, it is easy to see that the last equation would be 
modified as 



de 



cos 261 



(cos 61 -t 77)1/ 



jj^H{cos0 + 7^). 



(for |(H-l)|>eV4) 



(10) 



This function is shown in Figure |21 In fact, equation (|HJ|, 
together with ^ , is valid to third order in e for nearly all 
values of This is because, for small e, the appropriate 
limits of integration of equation |2Jl are approximately 
equal for an elliptical and circular source. The small 
changes to the integration limits introduce O(e^) correc- 
tions to a quantity that already includes O(e^) terms; the 
accuracy of this statement can also be verified directly in 
the numerical results displayed in Figure 13 below. Since 
changes to the integration limits are small, the difference 
in magnification between an elliptical and circular source 
can be calculated simply by taking the difference between 
the radius of an ellipse and circle, Ar = (r — p), times the 
area element rdrd^, weighted by the local magnification, 
and integrated over all angles such that cos 6 + > 0. 
However, this approximation breaks down for I77I ~ 1, 
where the limb of the source is near the caustic. First 
consider sources just exiting the caustic with rj — 1. 
For < (p < 45°, the limb of a circular source will exit 
the caustic first, whereas part of the limb of an elliptical 
source will still be inside the caustic. The situation is 
reversed for 45° < 4> < 90°, when the elliptical source 
will exit the caustic first. This effect is not accounted 
for in equation |(2Jl. This formalism also breaks down 
near 77 = 1 for similar reasons. As a result, equation |(2Jl 
cannot be used when |(|?7| — 1)| < e^/4. 
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Fig. 4. — Contours of constant RMS fractional deviation of the 
light curve of a uniform, elliptical source from an equal-area uni- 
form circular source, as a function of the ellipticity e and position 
angle </> of the source. We have assumed that 20% of the light is 
contributed by sources unrelated to the fold caustic images. The 
solid contours are exact, the dotted contours are the approxima- 
tion 0.3e^ cos 2(/). The heavy solid curve shows the locus of (e,0) 
where the deviation due to ellipticity vanishes. Contours are at 
equally-spaced logarithmic intervals of 0.5 dex. 

Since there are no terms containing odd powers of e, 
corrections to equation © are of order . Therefore, the 
approximation is very accurate for rather high values of 
the ellipticity. We find that it predicts the magnification 
to better than 1% for e < 0.3 and \{\'q\ - 1)| > eV4. 
Figure 13 compares the difference in magnification given 
by equations © and (0) , with the exact calculation using 
equation ^ for three different values of the ellipticity, 
e = 0.5,0.3, and 0.05. 

So far, we have parameterized the magnification in 
terms of ry, the perpendicular angular separation of the 
center of the source from the caustic. To convert to 
magnification as a function of time i, we can replace 
r/ — !■ (t — to) //St, where is the time when the cen- 
ter of the source crosses the caustic, and we have defined 
At = pt E CSC -f. Here ts is the time is takes the source 
to cross an angle 9e, and 7 is the angle the trajectory of 
the center of the source makes with respect to the caus- 
tic. One cannot constrain and 7 separately from a 
measurement of At using observations of a single linear 
fold caustic crossing, since is degenerate with 7. How- 
ever, with two crossings, which will generally have two 
different values of 7, both parameters can be measured. 

In terms of detectability, the quantity of interest is 
root-mean-square (RMS) fractional deviation a between 
an elliptical and circular source. Under the approxima- 
tion equation (jS)), this is simply 

cr(e,(/)) = Ae^cos2(/), (11) 

where A = (rya - mV^ ^ dT]E{r])/{Go{T]) + 0.2) is the 
RMS of the normalized fractional deviation due to el- 
lipticity over the range (771,772)- The total magnification 
of images not associated with the caustic is generally of 
order unity. For typical source sizes of p ^ 0.05, the 
mean total magnification of the fold images during the 



crossing (|ry| < 1) is ~ p^^^^ ^ 5. We have therefore as- 
sumed a 20% contribution from images unrelated to the 
fold caustic. For — 1 < 77 < 1.5, A ~ 0.3. We note that 
this value of A depends not only on the adopted value 
of the contribution from images unrelated to the source, 
but also on the time- interval (771, 772) over which the RMS 
is evaluated. Changing this interval will increase or de- 
crease the value of A, and therefore change the evaluated 
RMS fractional deviation. We have chosen the interval 
— 1 < rj < 1.5, because this essentially the largest inter- 
val where the deviation from ellipticity is significant (see 
Figure EJ), and therefore this is the interval of interest in 
terms of signal-to-noise. 

Figure 21 shows the RMS deviation in this range as a 
function of e and 0, for both the approximation in equa- 
tion (|ll(l . and the exact calculation. It is clear the the 
approximation is excellent for e < 0.3. It is interesting 
to note that the deviation due to ellipticity vanishes for 
some value of </), for all values of e. For e ^ 1, this occurs 
at (f) = 45°, whereas this null shifts to larger values of 
as e — > 1. Interestingly, for large ellipticities, the light 
curve from a circular source is degenerate with that from 
an ellipse with its major axis aligned nearly parallel to 
the caustic. 

It is clear from the form of equation (jS]) that one can- 
not measure the ellipticity of a source from one caustic 
crossing alone; rather one measures only the quantity 

cos 20. If the source crosses multiple caustics, each of 
which will generally have a different value of </>, then it 
will be possible to measure e and (f> separately. 

3. APPLICATIONS 
3.1. Stellar Oblateness and Limb Darkening 

Microlensing has been proposed and employed as a 
method of studying var ious topics in stellar astrophysics 
(see, e.g. lGouldll2001l for a review). The fold caustic- 
crossings created by binary-lens events toward the Galac- 
tic bulge and Magellanic clouds are exceptionally useful 
in this regard, as t hey are rel atively fre quent ^, and they 
can be anticipated ijJaroszvhs ki & Maall200lj) . 

Here we consider whether the oblateness of stars can 
be measured in caustic-crossing binary-lens events. The 
oblateness / of a star is defined as the fractional differ- 
ence between its equatorial and polar radii. This ca n be 
related to its rotation rate fe.g. lSeaeer fc Huj|l2002|) . 

where M, _R, and P are the mass, radius, and rotation 
period of the star, respectively. The ellipticity of the star 
is thene= [1- (/-f l)-2]i/2. 

We first estimate the distribut ion of ellipticities of gi- 
ant stars in the Galactic bulge. lOlecfi fToOO) has com- 
piled a catalog of miscellaneous periodic variable stars 
in the Galactic bulge found by the OGLE II microlens- 
ing survey, which includes the period, de-reddened color 
and magnitude of each star. Most of these stars are red 

^ Approximately 10% of all microlensing events are caustic- 
crossing binary- lens events ^Icock ct al. 2000; Udalski et aTJ 
120001^ . whereas less than a few percent of microlensing events are 
single-lens events in which the source is resolved. Therefore, the 
majority of events for which it is possible to resolve the source are 
binary events. 
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Fig. 5. — The distribution of projected ellipticities of 551 chro- 
mospherically active giant stars in the Galactic bulge, as inferred 
from their colors, magnitudes, and rotation periods measured by 
OGLE II l 01ech..l996,) . The top axis labels show the inferred el- 
lipticity, and the bottom axis labels shows the RMS light curve 
deviations that would occur (relative to a circular source) if the 
stars were to cross randomly oriented lensing caustics. 



(V — I > 1), and are in the giant branch or red clump. 
K^Iechl ((1990) argues that most of these red variables, 
which comprise 1% of all stars in this part of the color- 
magnitude diagram, are likely to be chromospherically- 
active (i.e. spotted) stars. Under the assumption that 
all are single stars, and that the measured period cor- 
responds to the rotation period of the star, we use this 
catalog to estimate the distribution of oblateness for such 
chromospherically-active stars in the bulge. Since these 
stars are likely bulge giants, they will have masses of 
M ~ Mq. We estimate their angular radii from their 
{V — I)o color and Iq magnitude, employing the color- 
surface brigh tness relat ion of van Belle (1999), as pre- 
sented in Albr ow et al.l (^00). We then determine their 
physical radii by assuming that they are at a distance 
of 8 kpc. We determine the oblateness and ellipticity of 
each star using equation H12() . 

The three-dimensional ellipticity of a source is not ob- 
servable in caustic-crossing events, rather only the ellip- 
ticity projected on the plane of the sky, e± = e(l— cos^ i). 
Here i is the angle between the line-of-sight and the rota- 
tion axis of the star. For a ensemble of stars, this angle 
will be distributed uniformly in cosi. The probability 
distribution of ej./e is thus, P{e±/e) — ^[l—{e±/e)]~^^'^. 

We note that there is a selection effect that introduces 
a bias in our determination of the distribution of pro- 
jected ellipticities, which we make no attempt to correct. 
Nearly pole-on stars will have smaller variability ampli- 
tudes, and therefore will be preferentially missed by the 
OGLE survey. Therefore, our assumption that the dis- 
tribution of cos i is uniform is not actually correct for 
this sample. We expect this bias to be relatively small, 
because systems have to be very close to pole-on to be 
missed (due to the excellent quality of the OGLE pho- 
tometry), and because this bias will largely be washed 
out by the intrinsic scatter in the amplitudes of the pho- 
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Fig. 6. — The Unes show the relation e = [— 0.5r sec 20]-'^/^, 
the ellipticity e required for a uniform, elliptical source to ap- 
proximately reproduce the light curve of a limb-darkened circular 
source, as a function of the limb-darkening parameter F, for po- 
sition angles of the elliptical source relative to the caustic normal 
of (/) = 0, 10°, 20°, 30°, and 40°. Sources with ^ = 45° produce 
no deviation due to ellipticity, and sources with </) < 45° produce 
deviations that are incompatible with limb-darkening. Values of 
e > 0.745 are not allowed for rotating stars, as they would be 
rotating faster than breakup speed. 



tometric modulation due to, e.g. variable numbers and 
distributions of spots. 

Figure |31 shows the resulting distribution of projected 
ellipticities for the 551 stars from the lOlechI (|l996) cat- 
alog. The median ellipticity is e = 0.07, with ~ 20% of 
stars have e > 0.2. We relate the projected ellipticity to 
the RMS light curve deviation a by integrating equation 
(|11|) over (j>, which yields cr ~ G.2e^. The median value 
is (T ~ 0.1%, with > 15% of stars expected to produce 
RMS deviations > 1%. 

Current microlensing follow-up efforts can regularly 
achieve single-exposure precisions of 1% on bright giant 
stars. Therefore, if the majority of giant stars in the 
bul ge have rotat ion periods that are similar to those in 
the lOlechI l|1996|) sample, then 10 - 20% of all caustic- 
crossing events should show deviations due to the stellar 
oblateness. Whether or not chromospherically quiescent 
stars should have similar rotation periods to active stars 
is not clear. Our analysis demonstrates that this ques- 
tion could be answered with a sample of ^ 30 precise, 
well-covered caustic-crossing events. 

One potential complication to the measurement of 
oblateness is limb-darkening. For a circular, limb- 
darkened source with radial surface brightness profile of 
the form 



I{r) = 1 - F <; 1 - - 



P 



1/2^ 



(13) 



the magnification of a limb-darkened source is 
(pVlbrow ct al. 1999), 

A{t^- F) = A{^^ 0) + [^^ ''^ nCini^) - Go(^)]. (14) 
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The form for the surface-br ightness profile given in equa- 
tion l(T^ was introduced in ijAlbrow et alJ l999'). It is has 
the same behavior as the usual linear limb-darkening pa- 
rameterization, but has the advantage that there is no 
net flux associated with the limb-darkening. The normal- 
ized limb-darkening function Gi/2{ri)~Go(ri) is shown in 
Figure 13 

It is clear that the form of the deviation due to the 
limb-darkening is qualitatively similar to the deviation 
due to oblateness for > 45°, such that a limb-darkening 
circular source with F ~ — 2e^ cos(2(/)) approximately re- 
produces the light curve of a uniform elliptical source 
with parameters e,4>- The factor of two is the ratio of 
the RMS deviations of E{ri) and Gi/2(?7) — Go{r]) for 
— 1 < ^ < 1- The value of e required to produce a given 
limb-darkening F is plotted in Figure|Sl for several values 
of (f>. Figure|21shows an example, for a limb-darkened cir- 
cular source with F = 0.5, and a uniform, elliptical source 
with e = 0.5 and (p — 90°. Other examples are shown 
in Figure 131 The degeneracy is only approximate, how- 
ever the two curves agree well except near the beginning 
and end of the caustic crossing. Note that, in general, 
the time when the source exits the caustic (near 77 ~ —1) 
will be different for circular and elliptical sources of equal 
area. 

Of course, stars will, in general, be both oblate and 
limb-darkened. We find numerically that for small el- 
lipticities, the light curve of an elliptical, limb-darkened 
source is well-approximated by the superposition of the 
effects of limb-darkening and ellipticity. In the opti- 
cal to near-infrared, limb-darkening parameters for typ- 
ical microlensing s tars are in the range F = 0.3 — 0.7 
ijFields et alJl200^ . and therefore limb-darkening effect 
will typically dominate over the effect of ellipticity. How- 
ever, given the extraordinary precision of recent limb- 
darkening measurements (^Abe et al. 2003; Fields et al. 
|2QQ3), it is not clear that the effect of oblateness can 
be ignored . In particular, the inconsistency claimed by 
iFields et~ l. ( 2003) between derived limb-darkening pa- 
rameters of the K3 III source of microlensing event EROS 
BLG-2000-5 and stellar atmospheric model model pre- 
dictions may be partly reconciled by allowing for stellar 
oblateness. Fortunately, in this case a spectrum is avail- 
able from which the v sin j of the star can be constrained. 
However, in the future, modelers should be aware of this 
possible contamination to limb-darkening measurements. 

How can oblateness be distinguished from limb- 
darkening? First, since the effects are not completely 
degenerate, this may be possible to distinguish between 
them simply from high-precision, single-color measure- 
ments. Second, the oblateness signal is expected to be 
achromatic, whereas limb-darkening depends strong on 
wavelength. Finally, one can use multiple caustic cross- 
ings to attempt to distinguish between limb-darkening 
and oblateness: a second caustic crossing can be antici- 
pated, and in general, it would occur at a different angle 
(j) from the first. 

3.2. Oblateness of Extrasolar Planets 

IGraff fc Gaudil (1200(1) and ILewis fc Ibatal 

demonstrated that close-in, giant planetary compan- 
ions to the source stars of binary microlensing events 
can be detected if the planet crosses a fold caustic, 
which magnifies the refiected starlight to detectable 



levels. Subseque ntly, l Ashton fc Lewis! ()200H) and 
iGaudi. Chang. fc H an (20031) explored the detectability 
of structures associated with planets using this method, 
such as the planetary phase, atmospheric features, satel- 
lites, and rings. These authors found that, while the 
planet itself may be detectable, as well as variations due 
to the planet phase and associated rings, all other fea- 
tures are likely to be undetectable with foreseeable tele- 
scopes. 

Can the oblateness of giant planets be d etected via 
this method? .Gaudi. Chang, fc Hani ()2003j) calculated 
the expected signal-to-noise ratio Qp of the primary 
planet signal for a typical event. They found Qp ^ 
150(L'/100m)2(a/0.05 AU)-^, where D is the telescope 
aperture and a is the semi-major axis of the planet orbit. 
From equation Ullfl . the signal-to- noise ratio Qe of the 
deviation due to ellipticity, averaged over all ip, can be 
related to Qp by 

Q.^0.2eU^30e^(^)'(^)-^ (15) 

Thus we can expect signal-to-noise ratios of less than 
unity for the deviation due to ellipticity, unless the ellip- 
ticity is quite large, e > 0.2. On the other hand, close-in 
giant planets with a < 0.2 AU are expected to be tidally 
locked to their parent star, and therefore their oblateness 
due to rotat ion are expected to be quite small, / < 0.3% 
or e < 10% (jSeager fc Huill2002|) . Planets whose rotation 
periods have not been synchronized with their orbital pe- 
riods may have large oblateness. For example Jupiter 
has an oblateness of / ~ 0.065, which corresponds to 
an eccentricity of e ~ 0.34. Unfortunately, the amount 
of reflected light, and therefore the signal-to-noise ratio, 
falls off as a~^, decreasing the signal-to-noise ratio of the 
planet signal by a factor of > 16 for tidally-unaffected 
planets with a > 0.2, and therefore rendering the dis- 
tant planets too faint for the oblateness signal to be de- 
tectable, even if they are rapidly rotating. 

3.3. Accretion Disks 

Numerous authors have considered the idea of us- 
ing microlensing of multiply-imaged quasars to resolve 
the surface brightness distribution of the central ac- 
cretion disk. Studies range from theoretical examina- 
tions of the feasibility and details of th e method itself 
(e.g.. Agol fc Krolik 1999; Grieger. Ka vser. fc Refsdal 
119881 IGrieger. Kavser. fc Schramm .IQgi]) . to detailed 
fitting of light curves of the lens Q2237-I-0305 (e.g., 
iGoicoechea et al."2003'; 'Kochanekl l2003i: iShalvapin et alJ 
1^02; Yonehara 2001). The majority of these studies con- 
sidered source models of face-on accretion disks. This 
assumption is generally not appropriate for standard 
geometrically-thin accretion disks: disks are more likely 
to be seen edge-on than face-on. Here we briefly con- 
sider the signature of inclined accretion disks on fold 
caustic-crossing light curves, and in the process uncover 
a degeneracy between the inclination and position an- 
gle for disks with self-similar intensity profiles. For sim- 
plicity, we consider only linear fold caustics. However, 
we note that this assumption may be inappropriate for 
typical quasar microlensing scenarios, which have rela- 
tively large sources, and optical depths to microlensing 
near unity ()Kochanekll2003l:lWvithe. Webster, fc Turneil 



GAUDI & HAIMAN 



7 








3 







2 


S-H 















1 










_ 1 1 1 1 1 1 1 i 1 1 1 1 [ 1 1 1 1 1 1 1 1 1 ^ 

R=inR^ - 




R„-25Rg _Z 


ii 1 1 1 1 rr^-uij 1 


R~iR, 

l~T— r-t— 1-4-^ — 1 1 ~ 




Time (At) 

Fig. 7. — (a) Normalized surface brightness profiles rivir) as 
a function of radius r/Rs in units of the Schwarzschild radius, for 
a standard thin accretion disk with inner radius 3Rs- Profiles are 
shown for three different values of the radius R,, where the tem- 
perature of the accretion disk is equal to hv/k: Ru = 10-Rs (solid), 
Rv ~ 25-Rs (dashed), and ~ 4i?s (dotted), (b) Magnification 
as a function of time in units of At = R^jve^ where Ve is the ef- 
fective source transverse velocity. The curves correspond to the 
surface brightness profiles in panel (a). We have assumed that the 
disk has inclination i = 30°, position angle relative to the caustic 
normal of = 15°, an Einstein ring radius of B^Dob = 100i?Si ^^'^ 
that the magnification outside the caustic is Aq = 3. 



120001) . Lensin g of an inclined disk by a point mass was 
considered by ijHevrovskv fc Locb 1997j). 

We consider a sta ndard, optically-thick, geo metrically- 
thin accretion disk l)Shakura fc Sunyaevll973(l . The disk 
radiates locally as a blackbody with temperature, 



-3/4 



\3RsJ 



-1/2 



1 1/4 



(16) 



where r refers to the radial distance of a (face-on) disk 
from the center, and we have assumed an inner disk edge 
of 3i?s- Here is the radius at which the local tempera- 
ture of the disk matches the wavelength of observations, 
i.e. where T[r) = = hv/k. Under the assumption 
of a standard thin-disk model, Ri, can be related to the 
mass of the black hole and the mass accretion rate. The 
surface brightness profile is then, 



J,(r)cx (e^"/^W-l)' 



(17) 



Note that we have ignored all relativistic and Doppler ef- 
fects. Figure shows the normalized surface brightness 
profile for the assumption of Ri, = 10 -Rs, which roughly 
corresponds to the value found bvlKochanck ( 2003) from 
an analysis of OGLE light curves of Q2237+0305 in 
the F-band, corresponding to rest-frame 2000A. We also 
show surface brightness profiles for R^, — 2'*/'^10i?s and 
Ri, = 2~''/^10i?s, i-e. for half and twice the frequency of 
observations, respectively. 

We calculate the light curves by first determining the 
image areas of a finite number of concentric elliptical an- 
nuli using the exact form for the magnification (Eq. (SJ . 
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Fig. 8. — Contours show the Icr confidence limits on the incli- 
nation i and position angle c/), for simulated observations of a thin 
accretion disk crossing a linear fold caustic, with input parame- 
ters i = 20° and (f) = 15°, shown as the cross. We have assumed 
100 measurements from —10 < A4 < 10, each with 2% error, and a 
background magnification of Aq = 3. The lines types correspond to 
light curves in Figure|7| and are for Ry = 10-Rs (solid), Ry ~ 25-Rs 
(dashed), and Rv ~ 4.Rg (dotted). 



The ellipticity of a disk with inclination angle i (where 
i = is face on) is e = (1 — cos^ i)^/^. We weight each 
annulus by the local surface brightness, which is simply 
given by equation (|17() . with r replaced by the semi-major 
axis of the elliptical annulus. The magnification is then 
the sum over all annuli of the surface-brightness weighted 
image area, divided by the surface-brightness weighted 
area of the source. We integrate out to r = 10-Ri,, beyond 
which we find that the contribution to the magnification 
is negligible. Figure [7| shows light curves for a standard 
thin accretion disk with inclination i — 30° and position 
angle relative to the caustic normal oi 4) — 15°. The 
units of time are At = Rn/ve, where Ve is the effective 
source transverse velocity. In addition, we have adopted 
parameters appropriate to Q2237-I-0305, namely an Ein- 
stein ring radius of 0eDos = 100-Rg, and a total mag- 
nification of images unrelated to the caustic of ^0 = 3 
(see, e.g. lKochanekll20C)3l) . For this system, and y-band 
observations, At ~ 60 days. 

In 321 we demonstrated analytically that the light 
curves of uniform, equal-area elliptical sources with e <ti 
1 are degenerate in e and (j) , such that one can only mea- 
sure the combination e^cos2(/). In fact, as could have 
been anticipated from Figure^) we find numerically that 
the degeneracy between e and (j) exists for arbitrary ellip- 
ticities, again provided that the sources are normalized to 
have equal area. In turn, this implies that the degeneracy 
between e and (f) exists for any source with surface bright- 
ness profile that is elliptical and scale-free. The scale-free 
requirement arises from the fact that the ellipses must 
have equal-area, which implies different values for the 
semi-major axes. Therefore the introduction of a scale 
in the surface brightness profile will break the degener- 
acy. For accretion disks, the degeneracy translates into 
a degeneracy between and R^. However, since the 
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profile I,y{r) of a standard thin disk fEa. ll7|l is not scale- 
free, this degeneracy is not perfect. The severity of the 
degeneracy is set by the radius of the inner edge of the 
disk relative to R^. If Rs/Ru ^ 1, the surface brightness 
profile will essentially be scale-free, and the light curves 
will be degenerate. However, if Rs/Ru ~ 1, the central 
hole in the surface brightness profile of the disk creates a 
feature in the light curve (see Figured). Larger inclina- 
tions will produce more pronounced features due to this 
central hole, thus allowing one to distinguish the light 
curve produced from different inclinations. We note that 
this degeneracy was likely present in the simulations of 
IAgo1 KrolikT l)1999l) . although it does not seem to have 
been explicitly recognized as such. 

To crudely quantify the severity of this degeneracy, we 
have simulated observations of a caustic-crossing event. 
We adopted the parameters of the light curves shown in 
Figure 13 namely i = 30°, cj) = 15°, 6'eDos = 100i?s, 
and Ao = 3. We assume 100 equally-spaced observations 
from —10 < At < 10, each yielding photometry in a 
single band with a fractional flux error of 2%. For the 
parameters of Q2237+0305, this corresponds to a sam- 
pling interval of '--^ 12 days. This sampling interval and 
photometric error are comp arable to that in the OGLE 
light curve of Q2237-h0305 fWoznia k et 311120001) . Fig- 
urelHlshows the resulting la confidence regions in the i, (p 
plane for the three different light curves. As expected, 
the degeneracy is most severe for i?^ = 25i?s, where incli- 
nations as large as 80° are consistent with the simulated 
light curve at the la level. This translates into a factor 
of (cos30°/cos80°)i/2 ^ 2 uncertainty in R^. Note that 
we have fixed all other free parameters at their input val- 
ues; allowing these parameters to vary would worsen the 
degeneracy. 

Since a measurement of Ri, can be used to constrain the 
mass of the black hole, it is important to reduce the un- 
certainty in i?^ if possible. This can be done by obtaining 
higher-precision observations, although we find that, for 
the same sampling rate, the precision of the individual 
measurements must be < 0.2% to reduce the uncertainty 
to below 20%. A more robust way of reducing this un- 
certainty is to observe at higher frequencies, where the 
ratio Rs/Rv is larger. For example, for R^, = 10i?s, the 
uncertainty in R^ is only ^ 10%. Higher frequencies 
are desirable for two additional reasons. First, since Ri, 
is smaller, the source is more compact, which results in 
a higher magnification as the source crosses the caustic, 
improving the signal-to-noise ratio. Second, Doppler and 
relativistic effects are larger for radii closer to the black 
hole. These effects impose asymmetries in the surface 
brightness of the disk, which further reduce degeneracies 
between the parameters, and may even enable a m ea- 
surement of the black hole spin llAgol & KrolikllTool . 

It is also important to note that we have assumed ob- 
servations of a single, linear caustic-crossing event. Ob- 
servations of multiple caustic-crossings, or observations 
of crossings of more complicated caustic geometries (i.e. 
parabolic folds, cusps, etc.), would likely remove this de- 
generacy. However, the additional complexity of such 
geometries may give rise to other complications and new 
degeneracies; addressing these issues is beyond the scope 
of the present paper. 

4. SUMMARY 



Fold caustics are germane to many applications of grav- 
itational microlensing. Previous studies have primarily 
focused on sources with circular symmetry. Here we have 
considered microlensing of elliptical sources by fold caus- 
tics. We considered only linear fold caustics, which are 
generally applicable when the size of the source is much 
smaller than the Einstein ring radius of the system, and 
prove a useful and analytic approximation to the mag- 
nification structure near real fold caustics. The total 
magnification of the two images produced near such a 
caustic is proportional to the square-root of the perpen- 
dicular distance to the caustic, and thus diverges as a 
point source approaches the caustic. The chief utility of 
fold caustics is that this divergence can be used to achieve 
high spatial resolution and large magnification of faint or 
otherwise unresolved sources. 

By convolving the magnification near a fold caustic 
with sources of arbitrary ellipticity, we computed the 
magnification of a source near a fold caustic of scale 
Uf as a function of its ellipticity e, position angle cf) of 
the major axis with respect to the caustic normal, area 
7rp^, and separation from the caustic rjp. We demon- 
strated that, for e ^ 1, the deviation due to ellipticity 
is simply {ur/ pY^^e^ cos2(f>E{r]), where E{r]) is simple 
one-parameter function that can be trivially evaluated 
numerically or expressed in terms of elliptic integrals. 
The next higher-order corrections to this expression are 
of order e"*, and thus it is accurate to < 1% for e < 0.3. 
From this expression, it is clear that the deviation due to 
ellipticity vanishes for (/) — 45°, and one can only mea- 
sure the combination of parameters cos 2(f> from obser- 
vations of a single fold caustic crossing. For e <C 1, the 
root-mean-square light curve deviation due to ellipticity 
in the range —1 < rj < 1.5 is 0.3e^ cos 20. Surprisingly, 
we found that the deviation due to ellipticity vanishes at 
some value of > 45°, for all values of e, with (p > 45° 
for e > 0.3 and increasing with increasing ellipticity. 

We considered three applications of microlensing of el- 
liptical sources near fold caustics. We first demonstrated 
that, if most of the microlensing source stars have rota- 
tion periods similar to chromospherically active, spotted 
stars in the bulge, then approximately ~ 15% should ex- 
hibit light curve deviations due to their oblateness (cor- 
responding to ellipticities e > 0.2) that are detectable 
with current observations. The deviation due to limb- 
darkening is qualitatively similar to that due to elliptic- 
ity, such that a uniform elliptical source with parame- 
ters e, (j) can approximately reproduce the light curve of 
a circular, limb-darkened source with F = — 2e^cos20. 
This may complicate the interpretation of ultra-precise 
measurements of limb-darkening with microlensing. We 
then considered the feasibility of detecting the oblateness 
of close-in, giant, planetary companions to the source- 
stars of caustic-crossing microlensing events. We found 
that planets which are sufficiently close to provide a de- 
tectable reflected light signal are also generally tidally 
locked, and therefore rotating too slowly to produce a 
substantial deviation due to oblateness. If tidal locking 
can be prevented to preserve an ellipticity of e > 0.2 even 
in the close-in planets, then the oblateness may be ob- 
servable, especially for fortuitous geometries (with caus- 
tic crossing occurring along the direction perpendicular 
to the major or minor axis). Finally, we considered the 
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resolution of the structure of a quasar accretion disks 
using microlensing by stars in the foreground lens of a 
multiply- imaged quasar. We showed that there exists 
a partial degeneracy between the disk inclination, scale 
length, and position angle for the orientation of the pro- 
jected ellipse relative to the caustic. For passbands pri- 
marily arising from the outer portion of the disk, this 
degeneracy can lead to a factor of two uncertainty in the 



determination the disk scale length, which translates di- 
rectly into an uncertainty in the black-hole mass, in the 
standard thin accretion disk model. Higher-frequency or 
higher-precision observations can reduce this uncertainty. 



Work for BSG was supported by a Menzel Fellowship 
from the Harvard College Observatory. 
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